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Abstract 

Motivated by recent experiments, we calculate both the average current and 
the current fluctuations for a metallic island which oscillates between two 
symmetric electrodes. Electrons can only tunnel on or off the island when 
it is close to one of the electrodes. Using a Master equation we investigate 
the accuracy of such an electron shuttle both analytically and numerically. It 
is shown that optimum operation is reached when the contact time is much 
larger than the RC-time. 

I. INTRODUCTION 

One of the most exciting prospects of Coulomb blockade is the realization of a current 
standard based on periodically transferring single electrons through the device JIIH]- In 
order to minimize higher order tunnelling processes which limit the accuracy, a chain of 
tunnelling junctions was used ||. Quite recently, a single electron shuttle mechanism which 
is based on mechanical oscillations between two electrodes was suggested by Gorelik et al. 
Q. The major advantage of a mechanically driven single electron transfer is that higher 
order processes are negligible. 

In the original suggestion of a mechanical single electron shuttle a metallic island is 
connected via organic molecules to two electrodes. One of the electrodes is positive and 
the other is negatively biased. For mechanically soft molecules self-excitation is expected 
to internally drive oscillations of the island which transfers charge at frequency /. Because 
of Coulomb blockade an integer number of electrons is loaded onto the grain close to one 
turning point, and the same number of electrons is unloaded close to the other. The number 
of electrons transferred per cycle is thus expected to be 



(N) = 2n max , n Tl 
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where V is the applied voltage and C the island capacitance relative to the electrodes. The 
average current is given by (I) = ef (N) . If the island is situated between the two electrodes 
there should be either n max or — n max electrons on it depending on which electrode it is coming 
from. 

A first step towards an experimental realization of a mechanical single electron shuttle 
was undertaken by Erbe et al. [H]. The authors designed a semiconductor device where a 
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mechanical clapper oscillates between two electrodes - i.e. the oscillation is now externally 
driven with the average number of electrons transferred in each cycle of order 7 ± 2 || . 
This experiment which works similar to a bell is still in the classical regime, i.e. Coulomb 
blockade is not relevant here. The aim in a future experimental setup with a metallic island 
on the tip of the clapper is to obtain an electron shuttle mechanism which is limited to single 
electrons by Coulomb repulsion. 

Since the tunnelling resistance R increases exponentially with the distance Ax between 
the island and the electrodes (R = exp(Ax/A)), tunnelling is suppressed exponentially 
when the island is in the middle between the two electrodes (both the oscillation amplitude 
and the distance between the two electrodes is much larger than A). For the same reason 
both cotunnelling and simultaneous tunnelling between the island and the electrodes can be 
neglected. In order for electrons to tunnel the amplitude of the oscillation has to be large 
compared with the typical tunnelling distance A. In the experiments by Erbe et al. the 
oscillation amplitude is of order 100 nm (compared to a few nanometers as suggested by 
Gorelik et al.) and thus much larger than A . 

In our present work we investigate the accuracy of a mechanically driven single electron 
shuttle using a Master equation approach |J. The Master equation is solved numerically and 
also analytically in a simplified form. In particular, we calculate both the average number of 
electrons transferred per period which in general is different from the simple result (|l|) and 

the root mean square fluctuations AN = yj ' (N 2 ) - (N) 2 . As one of our main results, it turns 
out that optimum operation with AN -C 1 is obtained for long contact times (to ^> RC) 
and sufficiently low temperatures (k^T <C %A. 



II. MASTER EQUATION 



In order to calculate the probability p(m, t) to find m additional electrons on the island 
at time t we use a Master equation. This technique has been successfully used to describe 
both quantum dots and metallic islands exhibiting Coulomb blockade @- flnj. In our case 
the island oscillates mechanically which leads to time dependent transition rates T(m, t) at 
the leads. Collecting gain and loss terms the Master equation reads 
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In a golden rule approach the tunnelling rates are of the form || 
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Since both R and the capacitance C in AE oc |r are time dependent now, with the domi- 
nant time dependence associated with the exponentially varying R(t), the tunnelling rates 
factorize into: 
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Here 
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is a strongly varying function of time, dominated by the exponential decrease of tunnelling 
when the island moves away from the electrode (£ max is the time where the island is at its 
closest point to the right electrode). The motion is taken as purely harmonic, i.e. x(t) = 
x max sm(ut). For the rates r^^(t) whose time dependence is determined by that of the 
capacitance only, we take the standard result for tunnelling rates when a voltage of — V/2 
is applied at the left electrode and a voltage of V/2 is applied at the right electrode: 



ri&(t) = - 
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where r = RR(t max )C R (t max ) and C%(t) = C R (t) + C L (t). For the left electrode the indices 
R and L have to be interchanged, V has to be replaced by —V and x(t) by —x(t) . 

For the numerical calculations the full Master equation was used. For simplicity, we 
assumed that the capacitance depends linearly on separation like Ci.it) oc xo+x 1 . The 
function g(t) is of the form 
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which is sharply peaked around t max . Since tunnelling through the left and right electrodes 
does not take place simultaneously (x max ^> A), the tunnelling at the left or the right 
electrodes can be treated separately. For the analytic calculations we consider only one 
electrode: F^^m, t) = , which is approached by an island with n electrons from the left, 
i.e. we are calculating the conditional probability p(m\n, t) to find m electrons on the island 
given that initially there where n electrons on the island. 

It has to be noted that when the island is close to the right electrode, C L (t) changes 
much slower than Cr(£) - independent of the precise model for the capacitance. As long as 
the applied voltage V is not too large (i.e. as long as VCi,(t max ) is of order e, 2e, . . . ) Cl(£) 
can be replaced by Cx(£max) • In order to be able to solve the master equations analytically, 
we replace the function <?r(£) by a step function gn(t) having the same height and area: 
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The width 2t corresponds to an effective contact time. As expected, the contact time 
decreases for increasing oscillation frequencies (x max ^> A , t max ^> to). 

The solutions p(m\n, t) of the Master equation can be used to calculate the probability 
Pa(tti) that m electrons are transferred: 
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Pa(w) =Y,Pi( k )p( k ~ n\k,t + t max ) , 



(9) 



where p\{k) is the initial probability to find k additional electrons on the island {J2kPi(k) = 
1). To calculate this probability, we consider a symmetric situation where the probability 
for the island to approach the right electrode with n electrons is the same as the probability 
for the island to approach the left electrode with — n electrons: 

Pi(n) =Pf(-n) , (10) 

where 

Pi ( m ) = Pi{n)p{m\n, t + t max ) . (11) 

n 

This leads to a set of linear equations for Pi(n) which in principle can be solved for 
arbitrary values of n max and t . However, for practical purposes, small values for n max 
should be used. From p&{n) we then calculate both (N) and AiV {(^N^ = J2n N k pA{N)). 



III. LOW TEMPERATURE LIMIT 



For very low temperatures k%T <C |r Coulomb blockade fixes the number of additional 
electrons n on the island between — n max and n max strictly. We will focus on the behaviour of 



the single electron shuttle for small integer values of n n 

With the approximation described in eq. (H) the Master equation at T 
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We solve eq. (pD for t max - t Q < t < t max + t Q 
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where At = t max - 
using the fact that 
following result: p( 
Together with eqs. 

fluctuations for T = . The results obtained from the analytical expression (|14]) and also 



Eq. ( |14D can easily be proved by induction 
lX , the normalization condition leads to the 

n m ax+lM m M)- 
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and (PD this can be used to calculate the current and the current 



from the numerical solution of the full master equation are shown in fig. [T]. Coulomb blockade 
is clearly visible: for low voltages no electrons are transferred. Above a critical voltage of 

= 1,5 = the analytic expressions for the 
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the shuttle starts to work. For n r , 
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average number of electrons transferred {{N)) and the mean square fluctuations ((AN) 2 ) 
are comparatively simple: 
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where a = exp ( — • As expected, in the limit of very short contact times to *C r (a — ► 1), 
no electrons are transferred whereas very long contact times to ^> t (a ^ 0) lead to the 
simple result ([!]). The fluctuations disappear for both short and long contact times. Their 
maximum (AiV « 0.686) is reached for £ ~ 0.92r as shown in fig. |2|. 



IV. FINITE TEMPERATURES 

For finite temperatures the number of additional electrons on the island is no longer as 
restricted as for T = 0. Here, the Master equation has the more complicated form: 

^-p(-n max +k\n,t) = ~rj hrrTs ^^^(-n^ + k-^n^t) - p (-n max + k\n,t) 
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where e = exp (~ 2c^Mfc^r ) anc ^ £ = £ 2<5 • For t ^> r and constant Cs (i.e. constant e and £) 



the stationary solution (Jy?(m|n, i) = 0, which is independent of n (pf (m) = p(m\n } t mSLX +to)) 
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is given by: 

p { (-n max + k)= V ~ 1 e k2 e, (18) 
with the normalization constant r\ = J2^L-oo ■ Together with eq. @ this leads to: 

t— k oo 

PA(2n max + k) = ^ T Y, efe^r . (19) 
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In reality however, C%(t) is time-dependent and thus the results flIB|) and flT9| ) are not exact. 
It turns out however, that if one replaces C^(t) by C%(t*) where t* is defined by 

£ /UJ g(t)dt = r (20) 

one obtains good agreement between numerical results obtained by the solution of the time 
dependent problem and the analytic result (|T8|). (Changes which happen on time-scales 
much smaller than - — i* i.e. effective time-scales much smaller than r, are too fast to be 
followed, cf. fig. 0.) 

Fig. |^ shows two steps of the symmetric resulting Coulomb staircase for several tem- 
peratures. In the middle of one of the steps of the Coulomb staircase (5 = 0) the average 
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number of electrons transferred per period agrees with the simple result (N) = 2n max even 
for k B T > ^r- ■ The root mean square fluctuations are approximately given by 



As expected, the fluctuations disappear exponentially for decreasing temperatures. 



Using a simple Master equation we have derived analytic expressions for both the average 
current and its root mean square fluctuations in a mechanically driven single electron shuttle. 
We find a good agreement between our analytical results and a numerical solution of the 
fully time-dependent Master equation. 

To minimize fluctuations the shuttle should be operated under conditions where the 
interaction time to is much greater than the i?C-time r (to ~ 10t)) and the temperature 
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k B T is much lower than the Coulomb energy 

We would like to thank R. Blick, A. Erbe and N. Chandra for useful discussions. Support 
by the SFB 348 of the Deutsche Forschungsgemeinschaft is greatfully acknowledged. 




k B T « ^ 
k B T » ^ 




(21) 



V. SUMMARY 
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FIGURES 
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FIG. 1. The average number of electrons transferred per period (left) and the root mean fluctu- 
ations (right) for T = 0. Coulomb blockade is clearly visible: up to a critical voltage ( yc<L ^ max ) _ l -j 
no electrons are transferred. The Coulomb staircase becomes symmetric for to » r . The agreement 
between analytical results (lines) and computer simulations (crosses) is excellent. 




12 3 4 

t /x 



FIG. 2. The average number of electrons transferred per period (solid line) and the root mean 
fluctuations (dashed line) for T = and vc ^ ma -^ = i. p or t Q <g t max no current flows, the 
optimum range is for to ^ *max where the fluctuations disappear. 
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FIG. 3. For to ^> t and T = (solid line) the Coulomb staircase is perfect. For increasing 
temperatures, the Coulomb-staircase is still symmetric, however it is washed out leading to an 
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ohmic behaviour for high temperatures. (The energy scale is given by E c = -^-) 
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